Spectral density and metal-insulator phase transition in Mott insulators within 

RDMFT 
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We present a method for calculating the spectrum of periodic solids within reduced density matrix 
functional theory. This method is validated by a detailed comparison of the angular momentum 
projected spectral density with that of well established many-body techniques, in all cases finding 
an excellent agreement. The physics behind the pressure induced insulator-metal phase transition in 
MnO is investigated. The driving mechanism of this transition is identified as increased crystal field 
splitting with pressure, resulting in a charge redistribution between the Mn e g and tig symmetry 
projected states. 

PACS numbers: 



Transition metal oxides (TMOs), the prototypical 
Mott insulators, are test-bed systems for new function- 
al within density functional theory (DFT) and many- 
body theories alike. Ground state spectra obtained from 
many-body theories are in good agreement with exper- 
iments. Moreover, the spectral density obtained using 
dynamical mean field theory (DMFT)[l| and the G W 
corr ected DFT [3 agree with each other even for subtle 
features such as symmetry and site projected spectral 
density. Single particle DFT spectra can also be made to 
agree with these many-body results by using two sepa- 
rate fitting parameters; the on-site Coulomb term U and 
the scissors shift A, where A is the difference between 
the experimental gap and the Kohn-Sham gap obtained 
using the LSDA+C/ functional^. 

Away from the ground-state TMOs show the rich 
physics of insulator-metal phase transitions. The classic 
Mott insulator, MnO, exhibits metalization under pres- 
sure. This phase transition is accompanied by a simul- 
taneous moment and volume collapses-Si- On the the- 
ory side, however, the physics of this phase transition 
is totally different for different methods; while DFT re- 
sults indicate that the increase in band width controls 
the phase transition [il, [l(J DMFT results, on the other 
hand, show that the main reason for metalization lies in 
the increased crystal field splitting [![. 

Recently, RDMFT has shown potential for corre ctly 
treating Mott insulators under ambient conditions! Ill 
RDMFT is an appealing alternative since it docs 
not require any system dependent parameters and thus 
is a truly ab-initio theory for treating strong correla- 
tions. However, it still remains to be seen how RDMFT 
performs away from ambient pressure conditions; can 
RDMFT capture the insulator-metal phase transition? 
What is the physics of this phase transition within 
RDMFT? In order to answer these questions one requires 
two things: (l)a magnetic extension of RDMFT and (2) 
information about photo-emission spectrum to shed light 
on the nature of the phase transition. The latter is a diffi- 



cult quantity to extract from RDMFT which, by its very 
nature, is a ground-state theory. 

In the present work we extend RDMFT to describe 
magnetic solids and further present a technique for cal- 
culating the photo-emission spectrum. We validate this 
technique by demonstrating an excellent agreement of the 
t2 g and e g resolved spectral density thus obtained, with 
the well established many-body methods like GW and 
DMFT. We further show that not only at ambient pres- 
sure but also away from it RDMFT correctly determines 
the spectra of Mott insulators and captures the physics 
of the insulator to metal phase transition. 

Within RDMFT, the one-body reduced density matrix 
(1-RDM) is the basic variable 
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7 (r,r')=A J d 3 r 2 ...d 3 r N ^(r,T2...r N )^*(r',r 2 ...r N ), 

(1) 

where denotes the many-body wavefunction and A is 
the total number of electrons. Diagonalization of 7 pro- 
duces a set of orthonormal Bloch functions, the so called 
natural orbitals[l3|, (j)^, and occupation numbers, n^. 
In the present work we have extended RDMFT to the 
truly non-collinear magnetic case by treating the natu- 
ral orbitals as two component Pauli-spinors, leading to 
the spectral representation: 7(1-, r') = J^ik nik< ^ ik ( r ) ® 
Vik( r ')j with the necessary and sufficient conditions for 
ensemble A-rcpresentability of 7 [HI requiring < < 
1 for all i and k, and X)ik n *k = A. 

In terms of 7, the total ground-state energy [3| of the 
interacting system is (atomic units are used throughout) 



m = - 



iti V J hm V 2 7(r, r') dV + J p(r)V cxt (r) d 3 r 
1 f p{v)p{v r 



■ d 3 r dV 



EM 



(2) 



where p(r) = tr CT7 (r, r), V C xt is a given external poten- 
tial, and -E X c we call the exchange-correlation (xc) energy 
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functional. In principle, Gilbert's [14( generalization of 
the Hohenbcrg-Kohn theorem to the 1-RDM guarantees 
the existence of a functional E[j] whose minimum, for 
fixed a V ext yields the exact 7 and the exact ground-state 
energy. In practice, however, the correlation energy is an 
unknown functional of 7 and needs to be approximated. 
While there are several known approximations for the xc 
energy functional, the most promising for extended sys- 
tems is the power functional[ll| where the xc energy is 



l7 (r,r')| 2 



with a indi- 



given by E xc [j] = -i f J d 3 r'd 3 i |r _ r>| 
eating the power in the operator sense. In view of the uni- 
versality of the functional E xc [7] , the value of a should, in 
principle, be system-independent. A few "optimum val- 
ues" of a have been suggested in the literature [ll, III 17 1. 
In the present work a is fixed to 0.656 for all materials 
studied. 

In order to devise a theoretical method to approxi- 
mately obtain the spectral density we start from the def- 
inition of the retarted Green's function written in the 
basis of the natural orbitals 



iG^,(t - t') = B(t t')(*Z\{a x (t),a{,(t')}\*»), (3) 

where A = {i, k} with the index i labeling the natural or- 
bitals for a given k, a, are the creation and annihilation 
operators associated with the complete set of natural or- 
bitals and \^q) is the neutral TV-electron ground-state. 
The spectral function A\\> (uj) can be written in terms of 
the Lehmann representation as: 

A xx ,(u) = -23G*,H = 2 7 r]T« , K|*f +1 ) (4) 
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S Af±1 |*f ±1 ) is satisfied by the ex- 



wherc H\^\ 

act (N ± l)-particle eigen-states of the Hamiltonian H. 
To deduce an approximate expression for the spectral 
function, we replace the complete set of cigen-functions 
±:L )} by the set of approximate eigen- functions ob- 
tained by adding/removing a single electron in a natural 
orbital to/from the exact correlated TV-particle ground- 
state: 
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(5) 



While these many-body states are clearly not a complete 
set, we do expect them to capture the dominant contri- 
butions to direct and inverse photo-emission. The set 
{Iff* 1 )} does not include (N ± l)-states where in ad- 
dition to adding/removing an electron, other electrons 
are excited from the ground-state, i.e. terms involving 
more than one creation/annihilation operator. Thus we 



do not expect to be able to reproduce plasmon satellites. 
Replacing in Eq. ((4]) the complete set of exact eigen- 
functions by this incomplete set of approximate eigen- 
functions (in Eq. ([5])) and using the fact that the natural 
orbitals diagonalize 7, i.e. (^g \a\,a^\^fQ) = S\^uq, we 
end up with the following approximation for the spectral 
function: 



A X y{uj) = 2nS X \> [n x 5{uj - ej + (1 - n x )5(uj + e+)] 



(6) 



with ef = Eg - E x 



We note that inspite of being 
approximate, the spectral function in Eq. (|6|) satisfies 
the exact sum rule, ^ Axx' {oj)duo = 1. Being, as 
function of u), a single 5 functon for each fixed A = (i, k), 
the spectral function in Eq. ((HJ) is reminiscent of a non- 
interacting mean-field type approximation. We empha- 
size that our approximation is by no means mean-field 
because the (N ± l)-states in Eq. ([5]) are correlated and 
may even be strongly correlated if, for example, I^qO 
represents the ground-state of a Mott insulator. 

It is a formidable task to determine quasi-particle life- 
times within the Lehmann representation: one needs to 
determine the position of an infinite number of peaks 
corresponding to an infinite number of N ± 1 eigen- 
states in Eq. (j4|) . This coalescence of peaks is described 
by an envelope function the width of which is propor- 
tional to the inverse quasi-particle life-time. In order 
to determine these life-times within the present formal- 
ism, one may use the more general form of the N ± 1 
states known from the extended Koopmans' theorem fljj ]. 

Ixf* 1 ) = Ea^J^ 1 )' where l^ 1 ) are given by Eq. 
([5|. This yields for the diagonal of the spectral func- 
tion, Axx(u) = 2ir[nxW x ~(u;) + (l-nx)W£(u>)] with 

M = Zj \4x\ 2 ^ ± [Eg - which has a 

finite width allowing to determine the quasi-particle life- 
times. Such investigation will be left to future. Within 
this article we wish to focus on the, so called, density 
of states, which is obtained by taking the trace of the 
spectral function (Ea-^ aa ) : 

DOS = 2ir^2 - e^) + (1 - n x )6(u + e+)] (7) 

A 

Now what remains is to calculate the excitation ener- 
gies ef = e£ = Eg - E^ 1 , where E^ 1 is the energy 
of the system with an electron, with specific momentum 
k, added/removed. While in experiments E^ 11 repre- 
sents the total energy of a macroscopic block of material, 
in the theoretical description E^ 1 is the total energy 
of a large but periodically repeated Born-von Karman 
(BvK) cell, where a constant charge background is added 
to keep the total (infinite) system charge neutral. The 
calculation of such total energies is computationally very 
demanding; requiring the number of k-points times the 
number of natural orbital (typically ^2500) constrained 
ground-state calculations. This is a formidable task and 
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hence we make another simplification which is not con- 
ceptual in nature but rather a numerical trick similar to 
the Slater transition state procedure [l9j]; we first intro- 
duce the total ground-state energy, E^ ±v , where a frac- 
tional number of particles, 77, has been added/removed at 
a given ik. We then assume that upon adding/removing 
charge to ik the only occupation number that changes 
significantly is the one that corresponds to the very same 
ik while all the other occupation numbers as well as the 
natural orbitals remain unchanged. With this simplifi- 
cation, following Slater [l9j], the e ± can be approximated 
as 



di] 



dE[{<f>},{n}} 



,7=1/2 



On 



(8) 



™sk=l/2 



This simplification can be easily numerically validated by 
plotting the total energy as a function of n^; we find a 
nearly linear behavior for all the materials studied here. 
This implies that the Slater-type evaluation of the total- 
energy difference in Eq. (JSj> is rather accurate. 

Following the above procedure the spectral density for 
the strongly correlated Mott insulators NiO, CoO, FeO 
and MnO is calculated using the full-potential linearized 
augmented plane wave (FP-LAPW) code Elk[20j], with 
practical details of the calculations following the scheme 
described in Ref . (111. 
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FIG. 1: (Color online) Density of states for the TMOs in pres- 
ence of AFM order. Site and angular momentum projected 
spectral densities are also presented for transition metal e g 
and t2 g states and Oxygen-p states. In addition XPS and BIS 
spectra (shifted up for clarity) are presented for comparison. 
Again, a = 0.656 for all materials. 

It is immediately apparent from Fig. [T] that RDMFT 
captures the essence of Mott-Hubbard physics: all the 
TMOs considered are insulating in nature. This fact 
was already noticed in the previous work [ll[ where the 
presence of a gap without any spin-order was deduced 
via a very different technique, namely the discontinuity 



in the chemical potential as a function of the particle 
number. A closer examination of the spectra for NiO, 
CoO and MnO reveals an excellent agreement between 
the RDMFT peaks and the corresponding XPS and BIS 
data. In fact, not only the peak positions, but also their 
relative weights are well reproduced. For the case of FeO, 
it must be recalled that Fe segregation, unavoidable in 
this compound, precludes the experimental realization of 
pure FeO samples. For this reason the only existing ex- 
perimental data are rather old and the presumably sub- 
stantially contaminated and broadened data present no 
distinct features that may be used for a clear comparison. 

The actual values of the insulating gaps that may 
be extracted from Fig. Q] are 4.5(4.3)eV, 2.6(2.8)eV, 
3.2(3.6)eV, 3.5eV for NiO, CoO, MnO and FeO respec- 
tively with the corresponding experimental gap given in 
parenthesis. The value of the local moments we find to be 
1.36(1.9)mb, 2.7(3.3)ms, 3.35(3.62)/i B and 3.38(4.7)p B 
for NiO, CoO, FeO and MnO respectively, again with 
the experimental values in parenthesis. There are two 
reasons for the smaller values of the magnetic moment 
within RDMFT compared to experiment. Firstly, the 
calculations are performed with the FP-LAPW method 
in which space is divided into spheres around the atoms, 
the so called muffin-tins, and the interstitial region. In 
the case of fully non-collinear magnetic calculations the 
magnetic moment per site is calculated by integrating 
the magnetization vector field inside the muffin-tin. This 
means a small part of the moment is lost in the inter- 
stitial. Secondly, the power functional induces a slight 
non-collinearity in the magnetization leading to yet more 
loss in the z-projected moment. 

In Fig. [1] we also present the site and angular momen- 
tum projected spectral density for the TMOs considered 
in this work. The electronic gap, as expected, always oc- 
curs between lower and upper Hubbard bands dominated 
by transition metal d-states. However, while for NiO one 
finds a significant component of oxygen-p states in the 
lower Hubbard band, for the other TMOs this hybridiza- 
tion between oxygen-p and TM-d states reduces, and is 
almost absent in the case of MnO, indicating that for this 
material the insulating state is driven mostly by Mott- 
Hubbard correlations. As a validation of our method for 
the calculation of spectra we may compare these features 
of the projected spectral density, and in particular the 
ordering in energy of the t2 g and e g states, with well es- 
tablished ab-initio many-body techniques such as DMFT 
and the GW method 0, Q. In all cases we find excellent 
agreement, signaling that the method we present here 
yields an accurate description of the detailed features of 
the spectral density. 

Theoretical methods used to study the ground-state of 
TMOs agree with each other as far as the spectral de nsity 
of TMOs at equilibrium volume is concerned (ll I2I [Tol 21 1, 



however, the actual values of the gap and the moment dif- 
fer depending upon the details of the calculations. How- 
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FIG. 2: Left panel: on-site magnetic moment (in /is) for 
MnO as a function of reduced volume (v/vo). Right panel: 
Mn d-band occupancy resolved into e g and t^g components. 
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FIG. 3: (Color online) t2 g (right panel) and e g (left panel) pro- 
jected spectral density for MnO as a function of energy (in eV) 
and reduced volume. Results are obtained using RDMFT. 
Spectral density for different reduced volumes are shifted ver- 
tically for clarity. 



ever, this agreement between various methods ends at 
ambient conditions, in that when attempting to study 
pressure induced insulator-metal phase transitions (IMT) 
the results vary wildly depending upon the method. In 
MnO experimental data point towards a first order IMT 
which is accompanied by a reduced volume (v/vq= 0.68 
to 0.63) and moment collapse (5 to 1 /ub)@|- DFT with 
LDA/GGA like functionals captures this collapse of vol- 
ume and magnetic moment and shows that the physics 



behind this phase transition is the simple widening of 
the Mn d-states due to increased itinerancy of the elec- 
trons at the reduced volume^ [T^]. On the other hand, 
DMFT shows that the phase transition occurs as a re- 
sult of the increased crystal field splitting with the width 
of the Mn d-states unchanged The picture obtained 
with correlated band theories is even more com plic ated: 
LDA+C/ shows a moment collapse but no IMT|2lJ, hy- 
brid functionals yield a phase transition to a semi-metalic 
and finally the self-interaction corrected DFT 
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state 

shows a transition to a metallic state with unusually large 



d-band width 21 



Given this wide spread of results, we study TMO phase 
transitions using RDMFT in the present work to shed 
light on this controversy. RDMFT is an ideal method 
for doing so as it is an ab-initio theory for treating 
strongly correlated systems not requiring any adjustable 
system dependent parameters. The results, obtained us- 
ing RDMFT, for the magnetic moment in MnO under ap- 
plied pressure is shown in Fig. [2] It is clear that the mag- 
netic moment collapses from 3.6/^s at optimal volume to 
0.54/is at a reduced volume. Further reduction of the 
volume does not change the moment. Within RDMFT 
we find a volume collapse of 11% which is higher than 
the experimental value of 6.6%. In order to investigate 
the prime reason behind this moment collapse we plot, 
in Fig. [2j the number of electrons in the Mn d-states as a 
function of volume. One notices a redistribution of elec- 
trons amongst the symmetry projected ti g and e g states. 
At the reduced volume of 0.8 the number of e g electrons 
(n eg ) starts to reduce finally leveling off at v/vo=0.711. 
This is accompanied by an increase in the ti g state charge 
(jHig)- This picture is fully concomitant with the previ- 
ous results obtained using dmftQ. In order to look 

e g states we 



at the detailed behavior of these ti g and 



have also plotted them as a function of volume in Fig. 
[3] The e ff -states move above the chemical potential as 
the volume is reduced, while the ^-states move below, 
finally ending up with a totally different spectral density 
in the metallic phase as compared to the Mott insulating 
phase. As a result of this rearrangement the crystal field 
splitting between the t^g and e g states increases. Despite 
RDMFT being a totally different approach as compared 
to DMFT, the symmetry projected spectral density as a 
function of volume looks very similar for the two meth- 
ods(see Fig. 3 of Ref. [l]). However, there exists a strik- 
ing difference between the two in that within RDMFT 
the Mn d-states widen at reduced volume as compared 
to those at v/vq = 1. This increase in the band width 
has its origin in a small shift in the spectral weight of the 
lower Hubbard band to lower energy. A close inspection 
of Fig. [3] shows that the change in crystal field splitting- 
is very large (~4.4eV) and certainly has a prominent role 
in the IMT in MnO. Band widening is a co-existing phe- 
nomenon which has a very small influence on the IMT. 
This band widening on the other hand drives the IMT if 
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LDA/GGA functional with DFT arc used. 

To conclude we have presented a method to calculate 
photo electron spectra within the framework of RDMFT. 
We have shown that the spectral information obtained in 
this way gives a detailed account of the strongly corre- 
lated nature of the TMOs, including the subtle interplay 
between Mott-Hubbard correlation and charge-transfer 
character in these materials. We validate this method 
not only by showing a good agreement of gross spectral 
features with experiments, but also by a detailed compar- 
ison of the angular momentum projected spectral density 
with that of well-established many-body techniques, in 
all cases finding an excellent agreement. We have fur- 
ther elucidated the physics behind the insulator to metal 
phase transition in MnO using RDMFT. For MnO the 
pressure induced phase transition is caused by the in- 
crease in crystal field splitting which in turn is the result 
of a re-distribution of charge amongst the states with t 2g 
and e g symmetry. The widening of the transition metal 
d-band is seen as a co-existing phenomenon, but is cer- 
tainly not the reason for the metalization. 
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